Relaxation Dynamics in Type-II Superconductors with Point-like and Correlated 

Disorder 
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We employ an elastic line model to investigate the steady-state properties and non-equilibrium 
relaxation kinetics of magnetic vortex lines in disordered type-II superconductors using Langevin 
molecular dynamics (LMD). We extract the dependence of the mean vortex line velocity and gyra- 
tion radius as well as the mean-square displacement in the steady state on the driving current, and 
measure the vortex density and height autocorrelations in the aging regime. We study samples with 
either randomly distributed point-like or columnar attractive pinning centers, which allows us to 
distinguish the complex relaxation features of interacting flux lines subject to extended vs. uncor- 
rected disorder. Additionally, we find that our new LMD findings match earlier Monte Carlo (MC) 
simulation data well, verifying that these two microscopically quite distinct simulation methods lead 
to macroscopically very similar results for non-equilibrium vortex matter. 

PACS numbers: 74.25.Uv, 74.40.Gh, 61.20.Lc 



I. INTRODUCTION 

Technical applications of type-II superconductors, es- 
pecially in high-field configurations, require effective pin- 
ning mechanisms to prevent flux creep and flow and to 
thereby avoid dissipative losses. A disordered system of 
vortex lines at a finite temperature, subject to pinning 
from randomly distributed point-like or correlated pin- 
ning sites, forms a remarkably complex system displaying 
a wealth of features. Naturally occurring weak point-like 
disorder already destroys the low-temperature Abrikosov 
lattice present in a clean system in favor of a new Bragg 
glass phase with quasi long-range positional order. The 
first-order melting transition^— is then replaced by a con- 
tinuous transition into a vortex glass phased—, in which 
translational order is completely lost. Hence, vortex mat- 
ter with weak, randomly placed, point-like disorder al- 
ready shows a very rich phase diagram^" 9 .. The intro- 
duction of correlated disorder, such as columnar pinning 
sites, results in a strongly pinned Bose glass phase with 
localized vortex lines and a diverging tilt modulus^—, 
which is accessible to analytical treatment via a mapping 
to the propagation of bosons in imaginary tim o 13 ' 14 . 

Since Struik's original investigation on physical aging 
in various materials- 5 , many glassy systems have been 
found to show physical agingi£. More recent studies 
confirm that many other systems show the character- 
istics of glass-like relaxation and agingil~— . Experi- 
mentally, Du et al. detected evidence of physical ag- 
ing in disordered vortex matter by demonstrating that 
the voltage response of a 2H-NbSe2 sample to a current 
pulse depends on the pulse duration^. Thin-film Monte 
Carlo relaxation and aging studies of a coarse-grained 
two-dimensional model were performed by Nicodemi and 
Jensen^Ir. 24 . Bustingorry, Cugliandolo and Dommgucz 
investigated the relaxation of vortex matter employ- 
ing Langevin molecular dynamics (LMD) for a three- 
dimensional line model, finding clear indications of phys- 
ical aging in two-time quantities such as the density- 



density autocorrelation function, the linear susceptibility 
and the mean-square displacemen t 25 ' 26 . 

In this paper, we report the results of a study com- 
paring the non-equilibrium relaxation kinetics of vortex 
lines in the presence of randomly distributed point-like 
disorder and correlated columnar disorder. The origin of 
point-like pinning sites can be either naturally occurring, 
or artificially introduced crystal defects such as oxygen 
vacancies. Similarly, correlated columnar disorder ap- 
pears either as line dislocations or in the form of mate- 
rial damage tracks stemming from high-energy ion irradi- 
ation. It is well-established experimentally that colum- 
nar disorder yields considerably enhanced pinning effi- 
ciency over uncorrelated point-like disorder—. Since lin- 
ear pinning centers are extended along one spatial dimen- 
sion, one also expects profound differences in the out-of- 
equilibrium relaxation of magnetic flux lines as compared 
to samples with randomly distributed point pins. Nat- 
urally, these differences can only be addressed in a fully 
three-dimensional model and numerical study. 

We perform extensive Langevin molecular dynamics 
simulations on a coarse-grained elastic line model of vor- 
tex matter. For non-equilibrium systems, it is important 
to compare and thus validate different microscopic real- 
izations and simulation methods, in order to ascertain 
that the resulting macroscopic features stem from phys- 
ical properties of the system and not from artifacts of 
the specific algorithm. We first present our LMD data 
on the steady-state vortex velocity and radius of gyration 
for driven flux lines subject to point-like or columnar dis- 
order. To validate our simulation code, we compare these 
results for attractive point pins with earlier findings from 
Monte Carlo (MC) simulations. We proceed to systemat- 
ically investigate the complex non-equilibrium relaxation 
behavior of a system of initially randomly distributed 
and perfectly straight vortex lines via various two-time 
observables. For point-like disorder, we again compare 
our novel LMD results with previously published MC 
data3£. The main focus of this work is the distinct relax- 
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ation behavior of flux lines in the presence of randomly 
distributed columnar and point pinning centers. 

Our paper is organized as follows: In the next Section 
we define and explain the elastic line model as well as 
our LMD algorithm and discuss the values of the differ- 
ent system parameters. We then introduce the quanti- 
ties that we are using in order to understand the non- 
equilibrium properties of interacting vortex lines in the 
presence of different types of attractive defects. Sec- 
tion HIT] is devoted to a discussion of the steady-state 
properties. We use this regime in order to validate the 
different algorithms used for the study of our system. 
Section llVl presents our numerical results. In a system- 
atic study we disentangle the different effects due to the 
line tension, the vortex- vortex interaction, the pinning to 
defect sites, and the finiteness of the system. We discuss 
our main finding on how the different types of pinning 
centers, point-like and extended columnar defects, affect 
the non-equilibrium relaxation process. Finally, we sum- 
marize our results in Section fVl 



II. MODEL AND SIMULATION PROTOCOL 

A. Effective Model Hamiltonian 

We consider in the following a system of N vortex lines 
in the London limit, where the penetration depth is much 
larger than the coherence length. In order to model 
the dynamics of the system we employ a fully three- 
dimensional elastic line descriptio n 12 ! 29 . The Hamilto- 
nian of this system is written as a functional of the in- 
plane vortex line trajectories rj(z), where z denotes the 
direction of the applied external magnetic field, and con- 
sists of three competing terms: the elastic line energy, 
the attractive external potential due to disordered pin- 
ning sites, and the repulsive vortex- vortex interactions: 
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The elastic line stiffness or local tilt modulus is given 
by h « r- 2 e ln(A Qb /£ab), where T" 1 = M ab /M c rep- 
resents the effective mass ratio or anisotropy parame- 
ter, whereas A a & and £ a b respectively denote the Lon- 
don penetration depth and coherence length in the ab 
crystallographic plane. The in-plane vortex- vortex inter- 
action is given by V{r) — 2e K a (r / X ab ) , with the zeroth- 
order modified Bessel function K (essentially a logarith- 
mic repulsion that is exponentially screened at the scale 
A). The Nd pinning sites are modeled by randomly dis- 
tributed smooth potential wells of the form 
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U D (r, z) = 4-P S(z - z a ) 



1 — tanh 5 



b 



where p > is the pinning potential strength, and r Q 
and z a indicate the in-plane and z position of pinning 
site a. Lengths are measured in units of the pinning 
potential width bo. Energies are measured in units of eobo 
with eo = (^>o/47rA Q ;,) 2 , and the magnetic flux quantum 
4>q = hc/2e. 



B. Langevin Molecular Dynamics 

We employ a LMD algorithm to simulate the vortex 
line dynamics. To this end, we discretize the system into 
layers along the z axis. The layer spacing corresponds 
to the crystal unit cell size Co along the crystallographic 
c-directio n 29 i 30 . Forces acting on the vortex line vertices 
can then be derived from the properly discretized version 
of the Hamiltonian ([lj . We proceed to numerically solve 
the (overdamped) Langevin equation 
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dri(t,z) 5H[Ti(t,z)] 
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fi(t,z). 
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with the Bardeen-Stephen viscous drag parameter rj. The 
fast, microscopic degrees of freedom of the surround- 
ing medium are captured by thermal stochastic forcing, 
modeled as uncorrelated Gaussian white noise fulfilling 
(fi(t, z)) = and the Einstein relation (fi(t, z)fj(s, z')) = 
Ar]k^,T5ij8{t — s)5(z — z'), which guarantees that the sys- 
tem relaxes to thermal equilibrium with a canonical prob- 
ability distribution cx e ~ H / k n T . The time integration is 
performed via simple discretization of Eq. <j2j> 3 1 . 



C. Material Parameters 

We chose our simulation parameters to closely match 
the material parameters of the ceramic high-Tp type-II 
superconducting compound YBa2Cu307 (YBCO). The 
material is highly anisotropic with an effective mass 
anisotropy ratio of T -1 = 1/5. We set the pinning center 
radius to bo — 35A and measure simulation distances in 
terms of this length. The in-plane London penetration 
depth and coherence length are A a b = 34&o ~ 1200 A and 
£ a b = 0.3&o ~ 10. 5A respectively. The vortex line en- 
ergy per unit length is e ~ 1.92 • 10~ 6 crg/ cm, hence 
the line tension energy scale becomes t\ ps 0.189eo The 
depth of the pinning center potential is set to p = 0.05eo, 
except when noted differently. To fix the intrinsic sim- 
ulation time scale, we set the Bardeen-Stephen viscous 
drag coefficient rj — (f> l /2TTp n c 2 t; o 2 b « 10~ 10 erg ■ s / cm 2 
to one (for the normal-state resistivity of YBCO near 
Tc-i Pn ~ 500,uf2cm, see table 1 in Ref. HI]), resulting in 
a basic time unit of t — rjbo/eo ~ 18 ps. 



D. Relaxation Simulation Protocol 

Throughout our study, the investigated systems con- 
tained N = 16 vortex lines with L = 640 number of 
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layers (except where noted differently). In the scenarios 
that include disorder, the number of pinning sites per 
layer is Np/L = 1116 which corresponds to a mean in- 
plane distance of 9&o between pinning sites. In the case 
of randomly arranged point defects, the pinning site po- 
sitions are chosen anew for each layer, whereas for colum- 
nar disorder, each layer repeats the pattern of the first 
layer's randomly chosen positions. The system size is set 
to (16/V3)A a fc x 8\ a b, with the aspect ratio chosen such 
that a clean system with interacting vortex lines repro- 
ducibly forms a hexagonal Abrikosov lattice configura- 
tion after equilibration. We employ periodic boundary 
conditions in the x and y directions, and free boundary 
conditions along the z axis. 

Our initial out-of-equilibrium condition consists of a 
system of perfectly straight vortex lines, placed at ran- 
dom locations throughout the computational domain. 
Since the vortex line elements do not yet fluctuate (i.e. 
their distance from the vortex line mean in-plane position 
is zero), the internal vortex line configuration effectively 
is at zero temperature, hence, the start of the simulation 
at t = is similar to an up-quench to a finite temperature 
T = 10 K. This is in contrast to their random spatial dis- 
tribution, which is equivalent to an infinite temperature. 
We then let the system relax towards equilibrium until 
the waiting time t — s (typically in the range of 2 to 
4096) is reached, when we take a snapshot of the system. 
We proceed to calculate various two-time quantities (see 
section III El below) at logarithmically-spaced time inter- 
vals with a simulation end time that is ten times larger 
than the waiting time. 



E. Measured Quantities 

In the steady state of our driven flux line system, we 
measure the mean vortex velocity v by extracting the ve- 
locities of each line element in the direction of the driving 
force Fd from the time stepping algorithm in LMD and 
average over all line elements in the system. In MC, we 
take the average displacement in the direction of the driv- 
ing force Fd over 30 MC steps and calculate the mean ve- 
locity. Using Faraday's law, we can relate the vortex line 
mean velocity to an induced electric field E = B x v/c, 
which translates to a voltage drop across the sample. 
Similarly, the driving force is related to an applied ex- 
ternal current via the Lorentz force Fd = |j x O B/£?|. 
Hence, a driving force vs. mean vortex velocity graph is 
equivalent to experimentally determined current-voltage 
(I-V) characteristics. 

To quantify thermal spatial fluctuations along the vor- 
tex lines, we compute the vortex line radius of gyration 
r g = yj ((r^ z — fi) 2 ), i.e. the root mean-square displace- 
ment from the lines' mean lateral positions. The angular 
brackets again indicate an average over line elements as 
well as noise and disorder realizations. This quantity is 
expected to show a maximum at driving forces just below 
the depinning transition. 



To further accurately capture the relaxation and aging 
dynamics of out-of-equilibrium disordered vortex line sys- 
tems, we measure two-time correlation quantities. Since 
we wish to compare the relaxation behavior in LMD with 
previously measured MC data, we utilize the same two- 
time observables as in Ref. [28l : the height-height autocor- 
relation function of the vortex lines, the two-time mean- 
square displacement and the density- density autocorrela- 
tion function. All these quantities depend on two times, 
labeled in the following as s and i, with s < t. 

The roughness or height-height autocorrelation func- 
tion of the vortex lines is defined by: 

C(t,s) = ((r M (t) -ri(t)) (r itZ (s) -Ti(s))) , (3) 

where r i z (t) are the in-plane coordinates of line i at layer 
z at time t, fi(i) is the mean position of line i, and the av- 
erages are taken over all line elements as well as noise and 
disorder realizations. This quantity contains information 
about local thermal fluctuations of vortex line elements 
around the flux line's mean lateral position. In the case 
of free, non-interacting vortices, it can be mapped to the 
height correlation of growing one-dimensional interfaces 
(see Sec. UVA")) . 

The two-time mean square displacement, defined as 

B(t,s) = ((r i , z (t)-r i , z (s)) 2 ) , (4) 

measures the average square distance between a vortex 
line element's position at time s and a subsequent time 
t. This quantity provides data on the time evolution 
of the global structure of the vortex line configuration, 
in addition to the same local information contained in 
C(t,s). 

Finally, the two-time density-density autocorrelation 
function C v (t, s) is an observable that is measured by sav- 
ing a snapshot of the positions of all vortex line elements 
in the system at time s, calculating the radial distance 
rt(t, s) each vortex line element traveled between times s 
and t and determining the number n c of vortex line ele- 
ments for which this distance is smaller than a prescribed 
cutoff distance ri(t, s) < r c . The density autocorrelation 
is then C v (t,s) = (n c /NL). Throughout our study, the 
cutoff distance is r c = 0.056o- The density autocorrela- 
tion also contains information on the formation or decay 
of global structures; thus we expect it to generally follow 
the behavior of B(t, s). 

F. Comparison of Microscopic Algorithms 

In order to simulate the dynamics of elastic lines in 
a disordered medium, an appropriate microscopic al- 
gorithm has to be chosen. A considerable amount of 
work has been done using Metropolis MC implemen- 
tations of the model described above and variations 
thereo f 29 ' 30 ' 33 . This is usually based on performing a 
random (but spatially truncated^) step and evaluating 
the ensuing change in the system's energy AE. The 
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Metropolis rule then accepts the step with a probabil- 
ity P(AE) = nan{l,exp(-AE/k B T)}. Gotcheva et al. 
investigated the differences between a Metropolis and a 
continuous-time MC algorithm for a system of flux lines 
on a discrete lattice and subject to varying temperature 
and driving forc o 34 ^ 35 . The continuous-time update rule 
preserved positional order, while the Metropolis rule led 
to a disordered moving state, questioning the validity 
of the Metropolis algorithm for studies of driven vortex 
matter in lattice simulations. More recent studies demon- 
strated that positional order was preserved in off-lattice 
Metropolis MC of driven vortex matter—. 

It is crucial to investigate and compare different micro- 
scopic implementations of algorithms such as Metropolis 
MC and LMD simulations in a non-equilibrium setting. 
The choice of algorithm might introduce spurious effects 
that cannot be predicted a priori. In order to separate 
actual physical effects of the studied elastic line model 
from these artifacts, we performed a careful numerical 
comparison of LMD with earlier MC studies. Fast mi- 
croscopic degrees of freedom are modeled by the ther- 
mal force term in LMD, see Eq. ([2]). In equilibrium the 
noise strength is set by the Einstein relation (fluctuation- 
dissipation theorem, FDT). In out-of-equilibrium situ- 
ations, there exists in general no FDT-equivalent that 
would uniquely determine the form and strength of the 
noise correlations. Since the large-scale and long-time 
characteristics of Langevin stochastic differential equa- 
tions can be drastically influenced by the noise correlator 
properties, it is necessary to validate results by compar- 
ing to other numerical methods 30 . In the first part of this 
article, we therefore compare results from a MC study for 
both steady-state and relaxation properties of the vortex 
line model to data generated by the LMD algorithm. It 
should be noted that a direct comparison of time scales 
is difficult since the length of a MC time step is a dy- 
namically generated quantity, whereas in LMD the time 
step duration is a function of the material parameters. 



III. STEADY-STATE PROPERTIES 

We first employ LMD simulations for interacting flux 
lines in the presence of point-like disorder with disorder 
potential strength p, subject to a driving force Fd stem- 
ming, via the Lorentz force, from an external current. 
Results for the steady-state velocity and gyration radius 
are gathered in Fig. [TJ The red dots in Fig. [TJ indi- 
cate MC-generated data, while the solid lines were pro- 
duced using LMD with different values of p. It is quite 
clear from Fig. Ufa) that a pinning potential strength 
p = 0.05e in MC corresponds top« 0.04e in LMD. 
In MC, vortex line elements test a region with a ra- 
dius of 0.25&o around their current position for possible 
jump targets. It is conceivable that the pin energy bar- 
rier appears a bit smoother in MC since its width falls 
into the same length scale, which leads to the observed 
renormalization of the pinning potential strength. The 
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FIG. 1: (Color online) Steady-state (a) velocity and (b) gy- 
ration radius of TV = 16 vortex lines with a length of L — 20 
elements as a function of the driving force Fd in the presence 
of point pins with varying disorder potential strength p. The 
red circles show data from MC simulations with p = 0.05eo, 
while the other graphs display LMD simulation data with 
p ranging from 0.03eo to 0.05eo in steps of 0.005eo- In (c) 
and (d), the force axis is rescaled with the disorder potential 
strength. The velocity curves in (c) cross at Fd/p « 0.19 while 
the gyration radius maxima in (d) align around Fd/p ~ 0.12, 
with a slight bias towards higher Fd/p for higher values of p. 



slightly higher maximum of the MC gyration radius data 
for p = 0.05eo over the corresponding LMD curve with 
p = 0.04eo in Fig. |TJ(b) supports this argument, since 
vortex line elements are most likely trapped at a pinning 
site until they escape via a single jump. It is much less 
probable for any line element to escape the pin via mul- 
tiple successive jumps. In LMD on the other hand, the 
thermal force is only an added component on top of the 
(in this case stronger) driving and elastic tension forces. 
Hence, LMD yields smooth escape trajectories out of a 
point pin's binding potential, and the vortex line is con- 
sequently less rough, resulting in a smaller radius of gy- 
ration. 

Figures [He) and (d) show that the depinning force 
scales roughly linear with the pinning potential strength 
p, as expected 1 ^. The vortex velocity curves cross at 
Fd/p ~ 0.19, while the gyration radii have their max- 
ima around Fd/p ~ 0.12, with a slight systematic shift 
to lower Fd/p for smaller pinning strength values p. It 
should be noted that a true continuous non-equilibrium 
depinning phase transition occurs only at T = (and 
in the thermodynamic limit). The gyration radius max- 
imum in the vicinity of the critical depinning force may 
consequently be understood as the thermally rounded 
remnant of this zero-temperature phase transition 3 ^. 

Figures [2^a) and (b) display the mean velocity v and 
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FIG. 2: (Color online) Steady-state (a) velocity and (b) gyra- 
tion radius of N = 16 vortex lines with L — 20 as a function 
of the driving force Fd in the presence of columnar disorder, 
with a disorder potential strength varying between p = 0.03eo 
and p = 0.05eo in steps of 0.005eo- As expected, the radius of 
gyration displays quite different behavior for columnar pins 
as compared to point pinning sites, compare Fig.fjjb). In (c) 
and (d), the force axis is rescaled with the disorder potential 
strength. The velocity curves in (c) cross at Fd/p ~ 1.9, while 
the gyration radius maxima in (d) align around Fd/p ~ 2. 



the radius of gyration r g as a function of the driving force 
Fd of vortices subject to randomly distributed columnar 
pinning sites. Correlated disorder is much more effec- 
tive at pinning flux lines than point disorder—. This is 
reflected in a critical depinning force that is about an 
order of magnitude higher for columnar defects than for 
uncorrelated point pins of the same strength p (per layer; 
results shown in Fig. [1]). In Figs. [2Jc) and (d) the driv- 
ing force is again scaled with the pin strength p. The 
vortex velocity curves cross at Fd/p ~ 1.9, which is in- 
deed a factor of 10 larger compared to point pins. In 
the presence of columnar defects, a single flux line may 
be in one of the following four configurations^ (i) un- 
pinned, located away from any pinning sites; (ii) trapped 
at a single columnar pin for its entire length; (iii) form- 
ing a vortex half-loop where the elastic line is trapped 
at a single columnar defect, with the exception of an un- 
pinned section that extends away from the defect line: 
Depending on the relative strengths of the driving force, 
thermal noise, and the pinning potential, the unpinned 
part may either expand or retract. This state represents 
a short-lived saddle-point configuration; (iv) forming sin- 
gle or double kinks by being simultaneously trapped at 
two adjacent pinning columns. This state is rather long- 
lived but will ultimately decay into either the unpinned 
or completely trapped state. 

The gyration radius indicates which configurations are 



typically assumed by the vortex lines. For Fd = 0, most 
vortex lines are fully trapped. The radius of gyration of 
a trapped line is restricted by the pinning potential ex- 
tension bo, hence we observe a marked reduction in the 
value of r g with columnar pins over free, unbound lines. 
With increasing but below-critical Fd, r g decreases since 
the mean position of trapped vortex lines shifts from the 
center of the pinning site, which further constrains fluc- 
tuations. In the vicinity of the critical depinning force, 
in the flux creep regime, r g rises sharply due to the for- 
mation of half-loops, single-, and double-kinks. Near the 
transition to free-flowing flux lines, r g develops a maxi- 
mum and gradually decreases for even higher Fd- In this 
state, vortex motion is slightly restricted by pinning cen- 
ters, but the flux lines move essentially unimpeded, and 
the radius of gyration approaches its unbound value. 



IV. RELAXATION PROCESSES 

In order to study the relaxation dynamics and possi- 
ble aging scaling of a system of vortex lines in various 
scenarios, we follow the procedures outlined in Ref. HH. 
As a test case for our simulation code, we first investi- 
gate a system of non-interacting lines without disorder, 
which can be mapped to the one-dimensional Edwards- 
Wilkinson (EW) interface growth model. We then pro- 
ceed to non-interacting lines in the presence of point- 
like disorder with two different pin potential strengths, 
where the short-time behavior is similar to the clean sys- 
tem while the long-time relaxation is modified by the 
attractive defects. To disentangle the effects of the mu- 
tual vortex repulsion from the disorder influence, we next 
study a system of interacting vortex lines without pinning 
sites. Subsequently, we present data on the full system 
of interacting vortex lines in the presence of point pin- 
ning centers. We then point out the differences in the 
relaxation kinetics in systems with point and correlated 
extended pins by investigating both non-interacting and 
interacting flux lines in the presence of columnar defects. 
Finally, we discuss finite-size effects due to short vortex 
lengths for both types of pinning sites. 

In each of the scenarios presented below, we first look 
at the relaxation of the single-time mean-square dis- 
placement B(t, 0), the time-dependent squared radius of 
gyration, r g (t), and the associated effective exponents 
P B (t) = d\nB{t,0)/d\nt and (3 h (t) = d\nr 2 g {t) / dhit. 
The mean-square displacement predominantly probes 
changes in the average positions of single vortices. Hence 
we use the average of /3b over an appropriate time in- 
terval as the aging exponent that we utilize to achieve 
(approximate) data collapse of the two-time global mean- 
square displacement B(t, s) and density autocorrelation 
function C v (t,s). The time-dependent radius of gyra- 
tion describes internal thermal vortex line fluctuations 
and enables us to compute an averaged scaling exponent 
0h that can be used for obtaining data collapse of the 
two-time height autocorrelation function C(t,s). 
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A. Free Non-interacting Vortex Lines 

The thermal fluctuations of the segment locations of 
free (non-interacting and not subject to disorder) di- 
rected elastic lines around their mean in-plane position 
can be mapped to the problem of a one-dimensional in- 
terface growing via random deposition. The continuous 
version of this growth model is described by the stochas- 
tic Edwards- Wilkinson (EW) equation^ 7 .. The tempo- 
ral evolution of the interface height relative to its mean 
height is governed by a diffusive term as well as a random 
noise term. Hence, it may be described mesoscopically 
by a linear Langevin equation, 



dh(z,t) d 2 h{z,t) 



dt 



dz 



■T){z,t) , 



(5) 



where v is the diffusive strength and rj(z,t) represents 
thermal white noise with zero mean and second moment 
(r](z, t)rj(z', s)) = 2TvS(z — z')S(t — s) that satisfies Ein- 
stein's relation. The temperature T enters through the 
noise strength. 

The linear nature of Eq. ([5]) makes it possible to 
arrive at analytical expressions for various two-time 
quantitie o 26 ' 38 ' 39 . In particular, the solution for the two- 
time height-height autocorrelation function in the corre- 
lated growth regime reads 



'58 



C(t,s) = C s 1/2 



1/2 



(6) 



Comparing with the general scaling form C(t,s) — 
s~ b fc{t/s) (with scaling function fc), this predicts the 
universal aging exponent b = 1/2 in the EW regim e 28 ' 38 . 
The mean square displacement follows a similar scaling 



form B(t, s) 



/b(V s )) wn ile the density autocorre- 



lation empirically scales as C v (t,s) = s fc„(t/s). 

The relaxation of the observables B(t, 0), r 2 {t) and the 
three two-time correlation functions B(t,s), C(t,s), and 
C v (t,s) in this scenario as observed in our LMD simu- 
lations is presented in Figs. |3] and |H respectively. For 
a very thin system with L = 10, we immediately see 
from the time evolution of the effective exponent /3b (t) 
[solid line with black circles in Fig. [3jb)] that the system 
starts to cross over into equilibrium where B{t, 0) ~ t, 
/3b — > 1 for early times t > 2 5 , and it is truly equi- 
librated after t > 2 10 . This is also visible in the un- 
sealed two-time quantities B(t,s) and C v (t,s) displayed 
in Figs. Hta,c): For s > 2 4 , the data of B(t,s) for dif- 
ferent waiting times s fall onto a single master curve, 
which indicates the recovery of time translation invari- 
ance. In contrast, the data for C v (t,s) collapse for all 
waiting times s. The difference in the onset of data col- 
lapse in these two quantities is caused by local thermal 
fluctuations contributing to the mean-square displace- 
ment, whereas short-scale variations are effectively av- 
eraged out due to the finite cutoff radius in the density 
autocorrelation 4 -. The effective gyration radius expo- 
nent /3/j(t), displayed in Fig. EJJd), decreases from the 




FIG. 3: (Color online) Relaxation behavior of (a) the vortex 
line mean-square displacement B(t,0), (c) the squared gyra- 
tion radius ri{t), and (b,d) the associated effective exponents 
/3b and j3h over time for free non-interacting flux lines with 
length L = 10 (• marker) and L = 2560 (► marker) averaged 
over at least 1000 realizations. 



start and eventually reaches zero around t ps 2 7 . For 
such a short flux line length, the crossover from the EW 
regime to the saturated (equilibrated) regime happens at 
very early times^ 8 ". The two-time height-height autocor- 
relation function C(t, s), plotted in Fig.@|b), shows data 
collapse for waiting times s > 2 7 , which reflects the time 
evolution of /3ft (t). 

The extended, bulk-like system with L = 2560 exhibits 
much slower relaxation and hence enables us to study 
the dynamical aging scaling regime. The effective expo- 
nents /3b (t) and j3h{t) [solid lines with black triangles in 
Figs. EIb,d)] show the remnants of a crossover for short 
times, while staying at a value of approximately 0.5 for 
t > 2 4 . Figures gjd-f) depict B{t,s), C(t,s) and C v {t,s) 
respectively as functions of t/s and scaled with appropri- 
ately chosen exponents of the waiting time. The data for 
all three quantities show data collapse for s > 2 4 with ex- 
ponent b — 0.5, indicating dynamical scaling and hence 
full aging in this time regime. The aging exponent of the 
two-time height-height autocorrelation function coincides 
with the predicted EW value from Eq. 

Except for the early-time crossover between the EW 
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FIG. 4: (Color online) Relaxation of two-time quantities in a 
system of non-interacting flux lines without disorder for a line 
length of (a-c) L — 10 and (d-f) L — 2560, averaged over 800 
different noise realizations. The sub-figures show (a,d) the 
mean-square displacement B(t,s), (b,e) the height autocorre- 
lation function C(t,s), and (c,f) the density autocorrelation 
function C v (t,s) as a function of (a-c) t — s and (d-f) t/s. 
Waiting times range from s = 2 to 2 12 = 4096. For L = 10, 
the system rapidly reaches equilibrium. The height autocor- 
relation function in (b) displays a short time window (s < 64) 
where it explicitly depends on s and time translation invari- 
ance is broken. For L = 2560 and s > 16, all quantities show 
aging and dynamical scaling with the EW exponent b = 0.5. 



and saturation regimes, which is not visible in Fig. 3 in 
Ref. HI, our findings produced via LMD simulations are 
in complete agreement with the MC data for both thin 
and extended "bulk" systems. 



B. Non-interacting Vortices with Point Disorder 

We now proceed to add point-like disorder with 
pininng potential strengths p = O.Oleo and p = 0.05eo 
to a system of non-interacting vortex lines, see Fig. [5] 
The equilibrium configuration at low temperatures con- 
stitutes an extremely dilute vortex glass. For the smaller 
defect strength of p = O.Oleo there exists an intermedi- 
ate time regime 2 3 < t < 2 7 during which both effective 
exponents and f3h(t) are fairly constant before de- 

veloping a maximum around t = 2 10 with a subsequent 
crossover into a frozen state with (3b, Ph where the 
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FIG. 5: (Color online) Relaxation behavior of (a) the flux line 
mean-square displacement B(t,0), (c) the squared gyration 
radius r^{t), and (b, d) the associated effective exponents (3b 
and (3h over time for non-interacting vortices subject to point 
pins with potential strength p = O.Oleo (• marker) and p = 
0.05eo (►), averaged over 1000 realizations. For p = O.Oleo, 
the dashed lines (below the curves) indicate the power laws 
with mean effective exponents (3b ~ 0.527 ± 0.003 in (a) and 
w 0.523±0.002 in (c) over the range 2 3 < t < 2 5 . Similarly 
for p — 0.05eo, the dashed lines (above the curves) show the 
power laws with mean effective exponents (3b ~ 0.725 ±0.041 
and #T w 0.716 ± 0.037 over the range 2 3 < t < 2 8 . 



vortex lines are firmly bound to the point defects [solid 
lines with black circles in Fig. HJb,c)]. The slight down- 
ward slope at early times is a remnant of the crossover 
from the random-noise into the EW regime of free flux 
lines. In fact, the exponent values at times t < 2 7 ap- 
proximately equal those of the disorder-free system, indi- 
cating that the time evolution of the mean lateral vortex 
position is essentially the same as for free lines prior to 
the disorder effects becoming noticeable. At later times, 
vortex movement starts to become affected by the at- 
tractive pinning sites, which temporarily accelerates the 
relaxation kinetics as the vortices are drawn into poten- 
tial wells, before flux line motion becomes at last frozen 
at the defects. 

FiguresUDja-c) show the resulting relaxation of the two- 
time mean square displacement B(t,s), the height-height 
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FIG. 6: (Color online) Relaxation of (a,d) the mean-square 
displacement, (b,e) the height autocorrelation function, and 
(c,f) the density autocorrelation in a system of non-interacting 
vortex lines of length L = 640, subject to randomly dis- 
tributed point pins with a potential depth of (a-c) p — O.Oleo 
and (d-f) p — 0.05eo; data averaged over 1000 realizations. 
Time translation invariance is broken throughout the simu- 
lation time window. For both pinning strengths, dynamical 
scaling for the mean-square displacement and the density au- 
tocorrelation approximately holds in an intermediate range 
of waiting times s, with the mean effective exponents inferred 
from Fig. [5] 



autocorrelation function C(t,s) and the density-density 
autocorrelation function C v (t, s) for the system with pin- 
ning strength p = O.Oleo- The global quantities B and 
C v are scaled using the aging exponent b — 0.543 taken 
from the average effective exponent /3b over the short- 
to intermediate-time region 2 3 < t < 2 7 , where [3b is 
roughly constant. These two-time autocorrelation func- 
tions yield approximate dynamical scaling for waiting 
times in this time regime, which further supports the in- 
terpretation that pinning sites are essentially irrelevant 
for the motion of the vortex line mean lateral positions 
during the early stages of the relaxation process. 

In the short-time regime t < 2 5 , the squared gyration 
radius rl{t) in Fig. [5jc,d) is also described by a power 
law with an approximate effective exponent f3h ss 0.523. 
Using this value as the aging exponent for the scaled 
two-time height-height autocorrelation function C(t, s) in 
Fig. in^b) reveals approximate dynamical scaling in this 
time window. However, stronger deviations from the free- 



line behavior are observed for C(t,s) than for the other 
quantities. For larger times t > 2 6 correlations become 
increasingly longer-lived, since vortex line elements be- 
come trapped at pinning sites. Hence, the influence of 
weak point defects is observed mainly in the fluctuations 
of flux line elements, whereas the movements of their 
mean lateral positions are hardly modified^. 

For a larger pin strength p = 0.05eo, the effective 
exponent maxima in Figs. [5jb,d) develop earlier, near 
t = 2 7 , and there appears no region with approximately 
constant exponents. For t > 2 11 , the system crosses over 
into a regime where the vortex line configuration appears 
to become frozen at the pins. This is reflected also in 
the two-time quantities, where we only see approximate 
data collapse for B and C v (t, s) for waiting times s < 2 6 
[Figs. [S)Jd,f)]. The two-time height autocorrelation func- 
tion C (t, s) in Fig. [6je) yields interesting non-monotonic 
behavior for s > 2 3 , where correlations actually increase 
again after developing a minimum. 

Elastic manifolds subject to disorder can be charac- 
terized by means of the roughness exponent x which is 
defined via the height-height correlation function along 
the manifold dimensions 43 (here, the contour length of 
the directed lines) 

C{z-z , ) = ([t(z)-t{z')]*)~\z-z'\ 3 *. (7) 

In the case of a dilute flux line system, for which mu- 
tual interactions may be neglected, and free of disor- 
der, thermal fluctuations lead to the EW roughness ex- 
ponent x = 0.5, in agreement with our numerical ob- 
servations. In the dilute vortex glass phase with point- 
like disorder, renormalization group analysis of manifolds 
subject to Gaussian disorder— ^ predicts a roughness ex- 
ponent x = 5/8. Our LMD simulations yield a distance- 
dependent effective roughness exponent in the range of 
0.5 < Xeff < 0.8 for a system of non-interacting vor- 
tex lines with point-like pinning sites. It should be em- 
phasized, though, that in our model the pinning sites 
are exclusively attractive. It turns out that the out-of- 
cquilibrium relaxation behavior is considerably different 
for directed lines subject to a mixture of attractive and 
repulsive pins: One then actually observes simple aging, 
albeit with non-universal scaling exponents that depend 
on temperature as well as pinning strengt h 28 ' 45 . 

C. Interacting Vortex Lines without Disorder 

To disentangle the effects of mutual flux line repul- 
sion from the influence of disordered point pinning sites, 
we next study the non-equilibrium relaxation of a clean 
system of interacting vortex lines. Fig. [7] shows the 
relaxation of the mean-square displacement, the radius 
of gyration, and the associated effective exponents in 
this scenario, whereas Fig. [8] displays the behavior of 
the three two-time autocorrelations. One may immedi- 
ately identify striking differences between non-interacting 
and mutually repelling vortex lines in the relaxation of 




FIG. 7: Relaxation behavior of (a) the flux line mean-square 
displacement B(t, 0), (c) the squared gyration radius r^(t), 
and (b, d) the associated effective exponents /3b and /3h over 
time for interacting vortices in a system without pinning 
centers, averaged over 5000 realizations. The dashed lines 
indicate power laws with the mean effective exponents (b) 
/3b « 1.57_± 0.08 (averaged over the time range 2 2 < t < 2 6 ) 
and (d) ^ « 0.50 ± 0.02 (averaged over 2 3 < t < 2 7 ). The 
inset in (c) shows the ratio of the x and y components of the 
radius of gyration for t > 2 12 . 



B(t,s) and its associated effective exponent com- 
pare Figs. |3ja,b) and Figs. [Zla,b). The initially large 
value of Pb in Fig. [TJb) can be traced to the rapid for- 
mation of long-range order due to repulsive vortex in- 
teractions. After the Abrikosov lattice has fully formed, 
flux lines perform confined random walks due to the effi- 
cient caging from neighboring vortices, which eventually 
leads to a low effective exponent /3 b ~ 0.2 for t > 2 12 . 
The data collapse in Fig. HJb) shows that an averaged 
Pb ~ 1-6 may serve as the effective aging exponent for 
the two-time mean square displacement B(t, s) for short 
waiting times s < 2 6 . 

The difference between interacting and non-interacting 
systems does not appear as drastic for the squared radius 
of gyration r^{t) as for B(t, 0), compare Figs. [3[c,d) and 
Figs.[7][c,d). After t > 2 10 , the Abrikosov lattice starts to 
form, and the repulsive forces due to neighboring vortex 
lines increasingly suppress transverse flux line wandering. 




FIG. 8: (Color online) Relaxation of (a,d) the mean-square 
displacement, (b,e) the height autocorrelation function, and 
(c,f) the density autocorrelation in a system of interacting 
flux lines without disorder and L — 640; data averaged over 
800 realizations. The left-hand panels show the unsealed log- 
log data, whereas data on the right-hand side are scaled by 
the waiting time s using the mean exponents from Fig. [7] 



For t > 2 12 , the gyration radius components along the x 
and y directions assume slightly different values owing to 
the anisotropic hexagonal vortex line arrangement, which 
in our rectangular system is always oriented along the x 
direction; see the inset in Fig. [TJ^c). For small waiting 
times s < 2 6 , the height autocorrelation data can be 
collapsed with the EW aging scaling exponent b = 0.5. 



D. Interacting Vortex Lines with Point Disorder 

We are now in a position to investigate the system 
of interacting vortex lines subject to attractive point- 
like disorder. As expected, the global time evolution, 
see Fig. [UJa,b), which is dictated by the mutual vor- 
tex repulsion, is hardly modified by the defects. The 
effective exponent /3b (t) displays essentially the same 
behavior as in the absence of pinning centers; compare 
Figs-inUb) and[7fb). Similarly, the aging exponents and 
the overall shapes of the mean-square displacement and 
the density-density autocorrelation in Fig. HOta.c) al- 
most match the simulation results without disorder. For 
short times t < 2 6 , the effective gyration radius expo- 
nent /3h(t) in Fig. [^d) is quite similar to f3h{t) in the 
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FIG. 9: Relaxation behavior of (a) the flux line mean-square 
displacement B(t, 0), (c) the squared gyration radius r^(t), 
and (b, d) the associated effective exponents /3b and fih over 
time for interacting vortices in a system with point-like dis- 
order of strength p — 0.05eo, averaged over 1000 realizations. 
The dashed line in (a) indicates a power law with the mean 
effective exponent /3b ~ 1.54 ± 0.08 (averaged over the time 
interval 2 2 < t < 2 7 ). 



non-interacting case with point pins, see Fig. [SJd). For 
longer times, repulsive forces alter the relaxation of Tg{t), 
which tends towards higher values, Fig. IUc). Hence, 
the global observables B(t,s) and C v (t,s) are influenced 
mainly through the presence or absence of vortex- vortex 
repulsion through the ensuing mutual caging. The local 
quantity C(t,s), on the other hand, better probes infor- 
mation on the disorder present in the sample. 

In MC simulations 28 , an interesting, non-monotonic 
behavior was revealed in the height-height autocorrela- 
tion function for a system of interacting vortex lines sub- 
ject to point-like disorder: The height autocorrelations 
displayed a pronounced maximum for small waiting times 
s and ln(t — s) « 5. Yet this feature is absent when the 
corresponding system is investigated with our LMD al- 
gorithm; see Fig. IToTb). In the present study, we have 
taken the vortex mass per unit length to be small and 
neglected the inertial term in the Langevin equation; see 
Sec. Ill Bl The algorithm used in Ref. [HI assumes a fi- 
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FIG. 10: (Color online) Relaxation of (a) the mean-square 
displacement, (b) the normalized height autocorrelation func- 
tion, and (c) the density autocorrelation in a system of in- 
teracting vortex lines with point pinning centers of strength 
p = 0.05eo and L — 640; data averaged over 800 realizations. 



nite displacement per MC step which generates an effec- 
tive mass. This in turn gives rise to oscillatory behavior 
at short times. This interpretation is indeed confirmed 
by LMD simulations that allow for a mass term in the 
Langevin equation, as depicted in Fig. 1111 With increas- 
ing vortex mass, C(t, s) shows damped oscillations on top 
of the monotonic time dependence of the overdamped, 
zero-mass case. 

We observe two-step relaxation behavior, typical e.g. 
for structural glasses, in the normalized two-time height- 
height autocorrelation C(t, s) when plotted as a function 
of the time difference t — s, shown in Fig. ITOTb). This 
behavior was also reported in the previous MC study 28 : 
A /3-relaxation regime where C(t,s) is hardly changing 
and displays time translation invariance precedes the ul- 
timate very slow decay. We attempted to fit a stretched 
exponential function to our long-time results, but could 
not achieve satisfactory agreement with our data. 



E. Non-interacting Vortices with Columnar Defects 

To compare the effects of uncorrelated point pins with 
those of extended, correlated defects on the flux line re- 
laxation kinetics, we now investigate a system of non- 
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FIG. 11: (Color online) The normalized two-time height au- 
tocorrelation function C(t, s) for different values of the vortex 
mass m in a system of interacting flux lines subject to ran- 
domly distributed point pins of strength p — 0.05eo, at wait- 
ing time s = 8, and averaged over 1000 simulation runs. The 
data for m > 0, for which C(t, s) displays non-monotonic and 
even oscillatory behavior, was generated using the Briinger- 
Brooks-Karplus integratoi^'^. 



interacting vortex lines in the presence of columnar pin- 
ning centers with a pinning potential strength p = 0.05eo- 
We start by first considering the case of non-interacting 
flux lines relaxing in the presence of columnar defects. 
Figure [T2l shows the relaxation curves for B(t,0) and 
r 2 g (t) with their associated effective exponents /3b (i) and 
f3h{t)- The time evolution of the mean-square displace- 
ment B(t,0) is slightly accelerated compared to the 
disorder- free case for times up to t = 2 9 ; see Fig. [3] This 
indicates that the initial trapping of vortices at linear 
pinning sites happens during this time regime. At later 
times, the effective exponent approaches the value of free 
flux lines /3 b = 0.5, due to unbound vortex line wander- 
ing. In fact, only ~ 4 to 5% of the vortices are pinned at 
the end of our simulation time window. 

As mentioned above, the average number of pinning 
sites per layer in our simulations is the same for point-like 
and columnar defects. Hence the combined lateral cross 
section along the z axis of all pinning sites is much larger 
for point than it is for columnar pins. Consequently, 
the probability for a vortex line to encounter a randomly 
placed columnar pin during what is essentially a ran- 
dom walk is much lower than in samples with randomly 
placed point defects. A flux line is therefore most likely to 
be initially captured by a single columnar defect (rather 
than multiple pinning sites, which would lead to vortex 
kink configurations) along a short length span, and sub- 
sequently becomes completely trapped at that pinning 
center; this yields rather small values for the final gyra- 
tion radius, see Figs. [T27c) and HJb). This is in stark 
contrast with samples containing uncorrelated point de- 
fects, where a single vortex line becomes captured by 
many pinning sites, which leads to a larger terminal ra- 
dius of gyration because the line is stretched in random 
directions between multiple pins; see Sec. IIVB1 Thus, 




FIG. 12: Relaxation behavior of (a) the flux line mean-square 
displacement B(t,0), (c) the squared gyration radius f"g(i), 
and (b, d) the associated effective exponents /3b and /3h over 
time for non-interacting vortices subject to columnar pinning 
centers with p — 0.05eo; data averaged over 10000 realiza- 
tions. The dashed lines indicate the power laws with the mean 
effective exponents ~$b~ » 0.672±0.037 and fa « 0.510±0.006 
over the time range 2 3 < t < 2 s . 



point-like pinning centers typically generate rough flux 
line configurations, whereas columnar defects straighten 
bound vortices. This fact can be expressed as an effec- 
tive upward renormalization of the elastic line stiffness, 
reflected macroscopically as a diverging tilt modulus for 
the entire vortex system in the pinned Bose glass phased. 

Comparing with the disorder-free system, the time evo- 
lution of the effective exponent of the radius of gyration 
(3h{t) in the columnar defect case is also different: The ex- 
ponent initially assumes a value 0.5 consistent with EW 
scaling. It begins to deviate from the EW value at t « 2 9 , 
while in the non-disordered case this decrease does not 
occur until t « 2 12 , see Fig. HHd) and Sec. llVA"! This 
is due to the fact that thermal line fluctuations inside a 
single columnar defect are confined to the pinning poten- 
tial well, which causes early saturation of the equivalent 
EW growth process. An appreciable number of kinks and 
double-kinks due to vortices trapped at multiple colum- 
nar pinning sites would presumably alter this relaxation 
behavior, but the occurrence probability of kinks is rather 
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FIG. 13: (Color online) Relaxation of (a) the mean-square 
displacement, (b) the height autocorrelation function, and 
(c) the density autocorrelation in a system of non-interacting 
flux lines of length L = 640 subject to randomly distributed 
columnar pins of strength p = 0.05eo; data averaged over 
1000 realizations. Scaling attempts with the averaged effec- 
tive exponents taken from Fig. [12] do not yield data collapse 
for B(t, s) and C v (t, s). One achieves better data collapse at 
early waiting times s for the height autocorrelations C(t, s). 




FIG. 14: Non-equilibrium relaxation of (a) the flux line mean- 
square displacement B(t, 0), (c) the squared gyration radius 
rg(t), and (b, d) the associated effective exponents /3b and /3h 
over time for interacting vortices in a system with columnar 
defects of strength p = 0.05eo , averaged over 1000 realizations. 
The dashed line in (a) shows a power law with the mean 
effective exponent /3b ~ 1.53 ± 0.11, averaged over the time 
interval 2 2 < t < 2 7 . 



small in our system, as explained above. At low temper- 
atures, our non-interacting flux line system in the pres- 
ence of correlated disorder forms a very dilute Bose glass, 
where the number of vortex lines is much less than the 
number of pinning sites. 

Figure [13] shows the relaxation of the three two-time 
autocorrelation functions for different waiting times s as 
a function of the ratio t/s. We obtain data collapse for 
the height autocorrelation function C {t, s) when scaling 
with the appropriately averaged effective aging exponent 
for early waiting times s < 2 9 and t/s < 2, which is 
consistent with the EW regime. The more global mean- 
square displacement B(t, s) and the density autocorre- 
lation C v (t, s) cannot similarly be scaled to obtain data 
collapse, because the effective exponent /3b (i) is never 
even approximately constant througout the entire ob- 
served time interval, as is evident in Fig. 112( b). 



F. Interacting Vortex Lines with Columnar Defects 

Next, we again turn on the repulsive vortex- vortex 
interactions. Similar to the disorder-free system and 
the samples with point pinning centers, see Sees. IIV CI 
and IIVD1 caging effects accelerate vortex motion, and 
the shape of the single-time mean-square displacement 
B(t, 0) and its associated effective exponent are hardly 
modified by the disorder; compare Fig. I14f a.b). At very 
large times t > 2 12 , B(t,0) becomes flatter and ap- 
proaches a plateau owing to the confinement of vortex 
lines by attractive defects. This effect is rather more pro- 
nounced for columnar pins than for uncorrelated point 
defects. In fact, the maximal values of B(t,0) and r^(i) 
are both smaller in the case of correlated disorder, which 
indicates tighter binding to the pinning sites. 

The radius of gyration also shows similar time evolu- 
tion trends as compared to samples with point-like dis- 
order; see Figs. [14Tc) and IHfc). The effects of attrac- 
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FIG. 15: (Color online) Relaxation of (a) the mean-square 
displacement, (b) the normalized height autocorrelation func- 
tion, and (c) the density autocorrelation in a system of inter- 
acting vortices of length L — 640 with columnar defects of 
strength p = 0.05eo; data averaged over 800 realizations. 



tive columnar pinning sites set in later than for point- 
like pins, owing to the aforementioned differences in the 
encounter probability for flux lines and pinning centers. 
The accelerated growth of the gyration radius for t > 2 9 
is due to the pinning at multiple sites and the subse- 
quent formation of kinks, here facilitated by the strong 
repulsive forces. The non-monotonic behavior of r^(t) at 
times t > 2 13 is caused by the decay of previously formed 
kinks. The density autocorrelation C v (t,s) in Fig. ITSf c) 
becomes flat for waiting times s > 2 s , which supports 
the interpretation that the vortex lines are essentially 
trapped by this time and the only remaining relaxation 
process is the decay of metastable kink configurations. 

Figure irST b) depicts the normalized height-height au- 
tocorrelation function of this system for different wait- 
ing times s. This function shows non-monotonic behav- 
ior, but in this situation it cannot stem from an effec- 
tive mass, as we checked. The appearance of the maxi- 
mum indicates a fundamental change in the lateral fluc- 
tuations. Although we do not yet fully understand this 
phenomenon, we tentatively relate this observation once 
again to the decay of kinks in the long-time limit. 
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FIG. 16: (Color online) Non-equilibrium relaxation of the flux 
line mean-square displacement B(t, 0) and the squared gy- 
ration radius r^{t) over time for non- interacting vortices in 
systems with point-like and columnar defects with identical 
pinning strengths p = 0.05eo for (a,b) L = 2 and (c,d) L = 10, 
averaged over 10000 realizations. 



G. Finite-Size Effects 

Effects due to the finite flux line length are best ana- 
lyzed in terms of the second crossover between the EW 
and saturation regimes for non-interacting vortices in 
the absence of pinning sites; see Sec. IIV Al (The first 
crossover between the random thermal noise and the 
EW regimes only depends on the EW diffusion constant, 
which here corresponds to the vortex line tension.) The 
time at which this second crossover occurs depends on 
the square of the line length, tc — (£/24&o) 2 27r/eo (in 
the limit of large L)^. Our choice of L = 640 for most 
of the simulation scenarios in this paper thus provides 
a sufficiently long time window tc ~ 2 12 to observe the 
competing effects of vortex interactions and pinning. 

Of particular interest is the value of L at which the 
relaxational difference between point-like and columnar 
pinning sites becomes apparent. Figure [TBI shows a com- 
parison plot of B(t,0) and r 2 g (t) for both columnar and 
point-like disorder for very short vortex lines with L = 2 
and L — 10. In a purely two-dimensional system with 
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L = 1, where the flux lines are reduced to point particles, 
there is obviously no difference between the two types of 
pinning sites. But already for L = 2 we observe differ- 
ences in the long-time evolution of B(t,0). The curve 
for columnar defects can be almost exactly reproduced 
by using point pins and halving the number of pinning 
sites. Hence the difference at L = 2 is largely due to the 
lower effective density of columnar disorder, as discussed 
in Sec. IIVE1 Yet this equivalence does not extend to the 
time evolution of the gyration radius r^(t), which dis- 
plays small but significant qualitative differences for the 
two defect types throughout the simulation time window. 

For L = 10, the long-time difference in B(t,0) be- 
tween columnar and point-like pinning sites can also be 
explained by the lower effective density of columnar pin- 
ning sites. But the deviations appearing at much shorter 
times reflect genuine physical distinctions in the pinning 
behavior and ensuing relaxation kinetics. 

V. CONCLUSION 

In this paper, we have investigated the differences in 
the non-equilibrium relaxation features between systems 
of magnetic flux lines in the presence of point-like and 
columnar disorder. Proceeding in a systematic way, and 
considering different limiting cases, allowed us to disen- 
tangle the distinct contributions originating from the at- 
tractive pinning centers, the repulsive mutual vortex in- 
teractions, and the line tension. 

We validated both our Langevin Molecular Dynamics 
simulation code and the Monte Carlo algorithm used in 
previous studies in a genuine out-of-equilibrium setting 
by comparing the steady-state vortex velocity and ra- 
dius of gyration as a function of an external driving force 
to results from Monte Carlo simulations. As discussed 
in Sec. Ill FL both these simulation methods need to be 
tested and validated when applied to non-equilibrium sit- 
uations. We found that the pinning potential strength in 
MC is slightly renormalized as compared to LMD due to 
the (inevitable) choice of a maximal MC step size. 

The introduction of columnar instead of point-like pin- 
ning sites dramatically changes the steady-state proper- 
ties. As expected, the critical depinning force is enhanced 
by approximately an order of magnitude. The radius 
of gyration is suppressed via vortex line confinement in 
columnar pinning sites for a driving force well below the 
critical depinning force. At the transition, partial depin- 
ning leads to the formation of half-loops and kinks in the 
vortex lines and thus to a sharp increase in the radius of 
gyration. At even higher driving forces, flux line motion 
is not influenced by pinning. 

We carefully studied the relaxation towards equilib- 
rium of a system of initially perfectly straight and 
randomly-placed vortex lines under various conditions by 
observing single- and two-time quantities to again com- 
pare the effects of uncorrelated point pins and correlated 



extended defects, and to further validate our LMD code 
against previously published MC results. We investi- 
gated the possibility of data collapse and, hence, a sim- 
ple aging scenario, by appropriately scaling our two-time 
quantities. We started with free, non-interacting vor- 
tex lines and showed that our results completely agreed 
with the MC data and the predictions from the Edwards- 
Wilkinson interface growth model. We then system- 
atically introduced attractive pinning centers and mu- 
tual repulsive vortex interactions. Caging effects due 
to vortex-vortex interactions lead to a considerable ac- 
celeration in the relaxation of global quantities, such as 
the single-time mean-square vortex displacement. A re- 
cent study revealed that the MC two-time height-height 
autocorrelation function for a system with interactions 
and point-like disorder displayed non-monotonic behav- 
ior (shown in Ref. l28h . Comparing with data obtained 
with an additional inertial term in our LMD algorithm, 
we argued that these oscillations stem from an effective 
mass generated by the introduction of a maximal MC 
step length. 

We demonstrated that the relaxation behavior of vor- 
tex lines depends crucially on the type of disorder. 
The vortex and Bose glass phases display complex non- 
universal relaxation features that are highly dependent 
on the material parameters. Once a deeper under- 
standing of the transient behavior has been established, 
detailed information contained in such time-dependent 
quantities could be used to characterize material proper- 
ties and specific samples. Point-like disorder binds vor- 
tices to many pinning sites at once, while columnar de- 
fects capture entire flux lines. When comparing these two 
defect types, one needs to take into account the difference 
in the effective pin density and thus the distinct probabil- 
ity of vortex line elements to become trapped. One may 
characterize the vortex glass phase through the roughness 
exponent \ °f the spatial height-height correlation func- 
tion along the strongly fluctuating flux line trajectory. 
In contrast, correlated linear defects effectively enhance 
the elastic line stiffness and hence straighten the trapped 
vortices in the Bose glass phase. 

We plan to expand our study to the transient proper- 
ties of driven vortex lines. The resulting relaxation then 
is towards a genuine non-equilibrium state in contrast to 
the relaxation towards equilibrium, which we presented 
in this paper. Other avenues for further investigations 
are the study of different and more realistic initial con- 
ditions, such as magnetic field or temperature quenches. 
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